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Ultrasound-driven oscillating micro-bubbles have been used as active actuators in microfluidic 
devices to perform manifold tasks such as mixing, sorting and manipulation of microparticles. A 
common configuration consists on side-bubbles, created by trapping air pockets in blind channels 
perpendicular to the main channel direction. This configuration consists of acoustically excited 
bubbles with a semi-cylindrical shape that generate significant streaming flow. Due to the geom¬ 
etry of the channels, such flows have been generally considered as quasi two-dimensional. Similar 
assumptions are often made in many other microfluidic systems based on flat micro-channels. How¬ 
ever, in this paper we show that microparticle trajectories actually present a much richer behavior, 
with particularly strong out-of-plane dynamics in regions close to the microbubble interface. Us¬ 
ing Astigmatism Particle Tracking Velocimetry, we reveal that the apparent planar streamlines are 
actually projections of a streamsurface with a pseudo-toroidal shape. We therefore show that acous¬ 
tic streaming cannot generally be assumed as a two-dimensional phenomenon in confined systems. 

The results have crucial consequences for most of the applications involving acoustic streaming as 
particle trapping, sorting and mixing. 


Among the vast amount of microfluidics systems that 
have been emerging since the late 90’s, those involving 
the use of acoustic fields to enhance flow control (i.e. 
acoustofluidics) have experienced a great development in 
recent years due to their practical simplicity, versatil¬ 
ity and potentiality for applications as particle separa¬ 
tion m, particle trapping [4] , acoustic manipulation of 
droplets using surface acoustic waves [3 IS] or acoustic 
centrifugation [7]. Direct piezoelectric actuation of mi- 
crochannel walls generates acoustic waves whose interac¬ 
tion with suspended particles through streaming flows or 
radiation forces has been used to concentrate, trap, or 
separate them through either acoustic streaming (inter¬ 
action with the fluid) or radiation forces (interaction with 
the particles) [SUTO]. However, other type of applications 
may require higher shear stresses or just a different range 
of frequencies than those achievable by standing acoustic 
waves. Instead of actuating directly on the solid bound¬ 
ary of the channel, an alternative is to actuate on the sur¬ 
face of a sessile microbubble that has been deliberately 
created in the system. The kHz-frequency oscillations 
of the bubble surface then generate a steady streaming 
flow of adjustable intensity. Such a technique is partic¬ 
ularly effective in soft-wall microfluidic systems such as 
PDMS (polydimethylsiloxane) devices, where waves in 
the substrate are strongly damped. Ultrasound-driven 
microbubble streaming flows has been studied for appli¬ 
cations to promote fluid transport m , manipulation and 
poration of vesicles [12], and more recently for applica¬ 
tions such as particle separation and trapping dsmi]. 

One of the simplest and most natural configurations to 
insert stable sessile microbubbles in PDMS microchan¬ 
nels is to create blind side channels. When the channel 
is filled with water, gas is retained in the pocket which 
then adopts a semicylindrical shape (Figure [^. Due to 
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Figure 1. Typical experimental set-up. Left: The PDMS 
channel and the piezoelectric actuator are bonded to a glass 
slide. Right: Top View of the semi-cylindrical bubble in the 
microchannel and visualization of the acoustic streaming flow 
through different experimental particle trajectories. 


the lithography-based micro fabrication, a 2D planar ge¬ 
ometry is imposed and a reasonable assumption is that 
the two-dimensional shape of the channel and its semi¬ 
cylindrical bubble shape damps the third component of 
the flow field, leaving a 2D flow profile. The assumption 
has been proven valid enough to study general aspects of 
the phenomenon as the typical frequency resonances m 
or the two-dimensional streaming flow by using asymp¬ 
totic models m- 

In this paper, we demonstrate through experimental 
measurements that particles do experience off-plane dy¬ 
namics in the vicinity of the micro-bubble, which reveals 
an unexpected degree of complexity in confined fluidic 
systems driven by acoustic streaming. In the next lines 
we will describe the dynamics of the complex flow, point 
out the consequences for applications, as well as deter¬ 
mine in which circumstances one can ignore the three- 
dimensional features of the flow. In order to inquire 
on the nature of the flow, a three-dimensional solution 
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Figure 2. Trajectories of four different particles. Top left: Projection of particles trajectories in the XY plane. Top right: 
Orthographic projection of the particle trajectories. Bottom: Projection of particles trajectories in the XZ plane 


of the Stokes equation showing the same trend can be 
constructed, which compares well with the experimental 
results. A more detailed description of the theoretical 
model will be shown in a companion paper m- 

Experimental Setup.- The microfluidic device is 
made of Polydimethilsiloxane (PDMS) by soft lithog¬ 
raphy: The PMDS components are blended at a 10:1 
ratio, degassed, and poured into a SU-8 mold to cross¬ 
link for 24 hours. The cured PMDS is then bonded to 
a glass slide via oxygen plasma treatment. Typical mi¬ 
crochannels (see Figure have a height D=100 pm, and 
a width H=1000 pm, while the blind channel holding the 
semi-cylindrical microbubble has a typical width W=80 
pm. The microbubble is driven by ultrasound actua¬ 
tion through a piezoelectric transducer (10 mm diameter, 
Physik Instrumente, Germany) glued to the glass slide, 
next to the PDMS channel (Figure [^. An amplified sinu¬ 
soidal signal is sent to the transducer through a function 
generator GW-Instek AFG-2125. All experiments shown 
in the paper have been made at a typical frequency of 
^ 20 kHz, appropriate for strong excitation of microbub¬ 
ble of the present size m- In order to acquire more data 
at sufficient resolution, the flow can be slowed down by 
choosing a lower voltage at the piezo transducer. 

3D partiele traeking.- The particle trajectories and 
velocities are measured using astigmatism particle track¬ 
ing velocimetry (APTV) [TSl [19]. APTV is a single¬ 
camera particle-tracking method in which an astigmatic 


aberration is introduced in the optical system by means 
of a cylindrical lens placed in front of the camera sensor. 
Gonsequently, an image of a spherical particle obtained 
in such a system shows a characteristic elliptical shape 
unequivocally related to its depth-position 2 :. The im¬ 
ages of the particles in the microfluidic chip are taken us¬ 
ing an inverted microscope Zeiss Axiovert in combination 
with a high-speed GGD camera PGO.dimax at recording 
speeds in the range from 500 to 1000 fps. The optical 
arrangement consists of a Zeiss Plan-Neofluar 20x/0.4 
microscope objective lens and a cylindrical lens with fo¬ 
cal length /cyi = 150 mm placed in front of the GGD 
sensor of the camera. Monodisperse fluorescent spheri¬ 
cal polystyrene particles with nominal diameters of 2 pm 
(Microparticles GmbH) were used for the experiments 
(Ppg = 1050 kg m“^). Illumination is provided by a con¬ 
tinuous diode-pumped laser with 2 W power at 532 nm 
wavelength. This configuration provides a measurement 
volume of 600 x 600 x 120 pm^ with an estimated uncer¬ 
tainty in the particle position of ±1 pm in the 2 :-direction 
and less than ±0.1 pm in the x- and ^-direction. The 
particle trajectories detected for this frequency vary in 
length and in time scale, and add up to more than 50 dif¬ 
ferent trajectories. More details about the experimental 
configuration and uncertainty estimation of the APTV 
system can be found in reference [19]. 

Results.- Figureshows measured trajectories corre¬ 
sponding to four selected particles. Each trajectory oc- 
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Figure 3. Detail of a typical particle trajectory. The green and blue paths depict respectively the inner and outer particle 
loops, in which the particle travels respectively away from and towards the mid plane z = 0. Top left: Projection of the particle 
in the ZX plane and in the XY plane (bottom left). Top right: Particle vertical distance from the system’s midplane z = 0, 
notice the jumps that the particle performs in the closest position to the bubble’s surface. The red discontinuous line shows 
the results from simulations. Bottom right: Radial distance r = y/x‘^ -\- y‘^ of the particle to the bubble’s symmetry in time 


cupies different volumes separated by two planes of sym¬ 
metry: the mid-plane ^ = 0 and the bubble’s sagittal 
plane ^ = 0. A third symmetry plane would be expected 
if the bubble were perfectly cylindrical. Most of the ob¬ 
served particles are confined in between the symmetry 
planes, but particles crossing the symmetry planes are 
also eventually observed. 

Figure [^illustrates a portion of a trajectory: particles 
follow an almost flat trajectory at large r but experience 
discrete “jumps” into a different plane every time the 
particle’s orbit reaches its perigee, i.e. when the parti¬ 
cle comes to its closest distance to the bubble’s surface 
Particle trajectories typically follow a periodic 
pattern similar to the one depicted in Figure The par¬ 
ticle starts close to the mid-plane 2 : = 0 with series of fast 
loops characterized by short periods, short apogees (far¬ 
thest position from the bubble’s surface Tmax) and large 
perigees (trajectory in green in Figure]^. When pass¬ 
ing through the perigee, the particle suffers small jumps 
away from the mid-plane. As the particle approaches 
the channel ceiling/floor, the orbit reduces the perigee, 
enlarges its apogee and consequently its period becomes 
longer (blue trajectory in Figure]^. Now, the particle 


jumps towards the mid-plane every time the it hits its 
perigee. As the particle approaches the mid-plane again, 
the orbit decreases its apogee, reduces its perigee, and 
will repeat the same pattern again. Since particles suffer 
strong velocity gradients, when choosing the experimen¬ 
tal settings one faces a dilemma: either to resolve the 
short time-scale particle dynamics or the long time-scale 
dynamics. Since our objective is to get a general view of 
the particle orbit, we chose the latter option. Therefore, 
note from the bottom plots in Figure that the selected 
experimental settings do not allow us to resolve experi¬ 
mentally the portion of the particle trajectory closest to 
the bubble and therefore the streamlines are not closed. 

Another remarkable symmetry becomes visible when 
observing the velocity field in detail in Figure as par¬ 
ticles approach the bubble at low angles 0, they are ad- 
vected into the bubble’s horizontal plane ^ = 0 at ^ < 60° 
and immediately after at 0 > 60° away from it. Once 
the particle leaves the bubble behind, there is a net dis¬ 
placement in its z—position whose direction and value 
depends strongly on the relative position to the bubble. 
Interestingly enough, the critical angle 0 ~ 60° at which 
the value of W changes sign corresponds also to the angle 
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at which Vr changes sign, and to the position at which 
Vo reaches its maximum value (Figure]^. It is interest¬ 
ing to note that in the horizontal plane 2 ; = 0, and 
Vr reach their maximum values, while Vz reaches a min¬ 
imum. In the data shown, the maximum Reynolds num¬ 
ber that a 2-micron-particle experiences is Rcp ^ 10“^, 
and therefore inertial effects do not play any role. How¬ 
ever, even though the flow is clearly viscous, stresses can 
be significant. For example, the maximum shear rates 
experienced at the top/bottom walls can be as high as 
500 1/s, while in the bulk they reach 1000 1/s. Note 
that mechanisms as shear-induced migration cannot be 
invoked to explain the particle dynamics due to the low 
particle seeding employed in our experiments. It was nor¬ 
mally assumed that the streamlines were two-dimensional 
and mostly closed [20] . However when observing the pla¬ 
nar projection, one can observe particle paths crossing 
streamlines. It was normally concluded that the particle 
was perturbed by the bubble surface oscillations, forcing 
the particle to change streamlines [14]. Although such 
a mechanism is valid and probably occurs as well, the 
present results show that apparent streamline crossings 
are also caused by the particle following a complex 3D 
streamsurface that has a toroidal shape (Figure]^. That 
explains why apparent streamline crossing occurs even for 
very small tracer particles. Interestingly enough, parti¬ 
cles are very rarely seen making a full and closed 3D tra¬ 
jectory. The most likely reason is that the time required 
to go through the full three-dimensional orbit (a full pe¬ 
riod) is much longer than our observation time. Still, 
streamline crossing cannot be totally neglected since the 
chances of a perturbation to occur while the particle trav¬ 
els are higher for longer periods. 

Discussion.- Secondary flows unavoidably appear in 
three-dimensional systems even at low Reynolds numbers 
when a streamwise component of vorticity is developed 
in the presence of a boundary layer. A classical exam¬ 
ple of secondary flow is the so-called tea leaf effect ED, 
which explains why particles concentrate in the bottom- 
center of a stirred tea cup. But they are also present 
when a pipe flow enters a bend or over an immersed ro¬ 
tating disc. In our particular case, in addition to this 
common phenomenon we must add the effect of the os¬ 
cillations along the axis of the cylindrical bubble, which 
cannot be neglected. These effects can be taken into ac¬ 
count by a modification of the two-dimensional bubble 
streaming solution, which is governed by a stream func¬ 
tion [El [17]. Due to the small Reynolds number 

associated with the streaming, it is appropriate to seek a 
linear superposition involving 2p2d and axial solutions of 
the Stokes equations satisfying no penetration conditions 
at the bubble and the walls of the microchannel, as well 
as the appropriate symmetries about x = 0 and z = 0 ob¬ 
served in the experiment. The lowest order axial solution 
(weakest radial decay) is axisymmetric and is governed 
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Figure 4. Particle Velocity distributions. Top: off plane ve¬ 
locity component Vz, note that they only become significant 
in the vicinity of the bubble’s surface. The arcs drawn in the 
top left Figure represents the section plotted in the top right 
plot. Middle: Azimuthal velocity component Vq. Bottom: 
Radial velocity Vr with maximum values at the sagittal plane 
and at approximately 0 = 30° (bottom right plots). 


by a Stokes streamfunction ^ax of the form 


'ipaxir, z) 
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( 1 ) 

Here, Ug is the characteristic scale of the streaming near 
the bubble surface (quadratic in driving voltage), a = 
IF/2 is the bubble radius, and represents a modified 
Bessel function of the second kind. The velocity field 
due to the superposition of 2D and axisymmetric modes 
is formally given by 


u(r, 0,z) = \/ X 


Caxi’axir, Z)0 


} 


( 2 ) 


where 0 < Cax 1 sets the relative strength of the axial 
flow and is used as a fit parameter between the theory and 
experiment. Figure shows that fluid trajectories under 
the linear superposition 0 qualitatively agree with ex- 
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perimentally measured particle trajectories over short 
ms) and long s) time scales. 

The results here shown may also have crucial conse¬ 
quences in stirring/mixing in microfluidic viscous flows. 
From Aref [22] and Ottino [23] we learned that, in order 
to enhance mixing in viscous flows, one must generate 
striations or tendrils thin enough to increase the value 
of the gradients and to allow diffusion to transfer mass 
in reasonable time scales. In this sense, 3D flows can 
achieve better mixing rates even in a passive way since 
they stretch and fold the liquid more efficiently, as has 
been shown by the Herringbone mixer m or the 3D ser¬ 
pentine [25] . The three-dimensional character of stream¬ 
ing flow from semi-cylindrical bubble contributes to the 
good mixing that can be already achieved in the liquid 
surrounding one single bubble. Mixing can be further 
enhanced by simply adjusting the bubble’s aspect ratio 
to enhance 3D effects, in addition to modulating its ac¬ 
tuating frequency or adding more bubbles in different 3D 
positions [15]. 

Conclusions and Outlook.- In this paper we have 
shown that the previously assumed 2D character of 
acoustic streaming by semi-cylindrical bubbles does not 
hold when particles approach the bubble’s vicinity. In¬ 
stead, 3D flow effects are present that have crucial conse¬ 
quences on the particle trajectories and therefore on the 
fluid flow surrounding the bubble. The mechanisms driv¬ 
ing such a flow must necessarily be connected with the 
presence of axial bubble oscillations and their necessary 
fulfillment of the boundary conditions. 
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